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Abstract 

Many morphogcnctic processes involve mechanical rearrangement of epithelial tissues that is driven by 
precisely regulated cytoskeletal forces and cell adhesion. The mechanical state of the cell and intercellular 
adhesion are not only the targets of regulation, but are themselves likely signals that coordinate develop- 
mental process. Yet, because it is difficult to directly measure mechanical stress in vivo on sub-cellular 
scale, little is understood about the role of mechanics of development. Here we present an alternative 
approach which takes advantage of the recent progress in live imaging of morphogenetic processes and 
uses computational analysis of high resolution images of epithelial tissues to infer relative magnitude 
of forces acting within and between cells. We model intracellular stress in terms of bulk pressure and 
interfacial tension, allowing these parameters to vary from cell to cell and from interface to interface. 
Assuming that epithelial cell layers are close to mechanical equilibrium, we use the observed geometry of 
the two dimensional cell array to infer interfacial tensions and intracellular pressures. Here we present 
the mathematical formulation of the proposed Mechanical Inverse method and apply it to the analysis of 
epithelial cell layers observed at the onset of ventral furrow formation in the Drosophila embryo and in 
the process of hair-cell determination in the avian cochlea. The analysis reveals mechanical anisotropy in 
the former process and mechanical heterogeneity, correlated with cell differentiation, in the latter process. 
The method opens a way for quantitative and detailed experimental tests of models of cell and tissue 
mechanics. 

Author Summary 

Mechanical forces play many important roles in cell biology and animal and plant development. In con- 
trast to inanimate matter, forces in living matter are generated by active and highly regulated processes 
within and between cells. The ability to directly measure forces and mechanical stress on the cellular 
scale within living tissues is critically important for understanding many morphogenetic processes but is a 
serious experimental challenge. The present work proposes an alternative approach based on the analysis 
of images that provide a visualization of cell boundaries in two dimensional epithelial tissues. The method 
uses the assumption of force balance within the epithelial layer to infer, on the basis of image-derived 
geometric data, the mechanical state of each cell. The proposed Mechanical Inverse method is illustrated 
on the analysis of two examples: the initial step of the gastrulation process in Drosophila embryo, and 
the process of neurogenesis in the developing avian cochlea. 

Introduction 

Genetics and biochemistry are central to all aspects of biological function. Physics is often less recognized 
but still important at many levels, everywhere from intramolecular to organismal scales. In particular, 
many important aspects of cell behavior depend directly and indirectly on its mechanical state defined 
by its interaction with neighboring cells and adhesion to the extracellular matrix [l}j3]. Cytoskeletal 
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mechanics and cell-cell adhesion determine geometric properties of cells [lj|4||6], as well as the dynamics 
of biological tissues |5|[7j-[T3"]. In plants, cells do not move, but the rigidity of cellulose membranes 



makes mechanical stress an obvious factor for cell division and proliferation 14 15 . It is known that 
animal cell proliferation also depends on substrate adhesion and the degree of cell confinement [2j|T6j- 



19 . It has also been demonstrated that (stem) cell differentiation is affected by substrate rigidity [20] 
More speculatively, mechanical feedback interactions have been conjectured a role in coordination of 
growth during development [TJ[9j[2Tj[22] . Mechanical transformation of epithelial tissue is of course itself 
central to many morphogenetic processes: gastrulation [7] and convergent extension [l], to name a few. 
Understanding how mechanical state changes in cells orchestrate morphological reorganization of tissues 
is an open problem and a subject of much current work [lj[7j[8] 

Our present understanding of the role of mechanics as one of the regulatory inputs into the cell is 
strongly impaired by the difficulty of quantitative characterization of the mechanical state (i.e. stress 
and deformation) of the cell. Among the available techniques are laser tweezers [23] and "traction force 
microscopy" [18,24 performed on cultured cells. UV laser ablation allows to mechanically perturb tissues 
|8| |25||26] on cellular scale with the time-lapse imaging of subsequent relaxation providing information on 
the mechanical state of the tissue. The ablation approach is widely used on live preps, for example, in 
the study of Drosophila embryonic development. Yet, this technique is definitely not a "non-destructive" 
one. 

On the other hand one of the major recent technical advances in developmental biology is the im- 
provement of live fluorescent imaging. These provide high quality time lapse movies of developmental 
processes, including interesting morphological transformations such as gastrulation and convergent ex- 
tension [7 p6p7] . The purpose of the present investigation is to explore what insight into the mechanical 
state of cells may be gleaned from a quantitative examination of high quality images of the type shown 
in Fig. 1A. Our goal is to use image analysis as a non-destructive approach to obtaining quantitative 
measures of stress in these systems. Similar strategy has been pursued by the recently proposed "Video 



Force Microscopy" (VFM) approach by Brodland et al 28 . Our approach will differ from VFM in its 



assumptions about mechanical state of tissue, in the parameterization of forces and in the way imaging 
data is utilized. 

Below we shall define a general model parameterizing the mechanical state of cells in two dimensional 
epithelial tissue and provide a computational method for inferring these parameters from the observed 
geometry of the cell array. We shall study the sensitivity of the proposed Mechanical Inverse (MI) 
method to errors in cell geometry and identify conditions under which robust inference is possible. We 
then illustrate the proposed MI method by applying it to the analysis of two different biological processes: 



cochlear neurogenesis 29 and ventral furrow formation 27 



Materials and Methods 

Model of epithelial tissue mechanics. 

Our approach is based on the assumption that epithelial monolayers are in an instantaneous mechanical 
equilibrium, characterized by a static balance of the forces acting at intercellular junctions. The second 
important assumption is that epithelial mechanics is dominated by the actomyosin cortices and inter- 
cellular Adherens Junctions 1| both localized at cell boundaries which form a visible two-dimensional 
web, as shown in Fig. la. Thus we assume that mechanical state of the cell can be described by effective 
tension at the interface and the hydrostatic pressure in cell interior. Yet, because cells can independently 
regulate their mechanical state, e.g. by modulating myosin activity or cell-cell adhesion, we allow for the 
possibility of each intercellular interface to have a different effective tension, T a b, and for each cell to have 
a different internal pressure P a (where a labels cells and ab labels the interface between cells a and b), 
as shown in Fig. ID. Mechanical equilibrium then corresponds to the condition that the forces acting on 
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each "vertex" fl (defined as a junction of three cells and therefore of three interfaces) add up to zero. 

Let fi and fj be the vertices belonging to the interface ab and let = rj — be the vector from 
vertex i to j. The force exerted by this interface on vertex i is 



F % = T ^ l + \( p a - Pb)4ep a , (1) 



where a labels vector components in the xy plane and ep a is the anti-symmetric tensor (e xy = —t yx = 1 
and e xx — e yy ~ 0). As shown in Fig. 2, this expression accurately represents the Young-Laplace balance 
between interfacial tension and the pressure differential across the interface P a — P b = n a bT a b, as long as 
the interfacial curvature n a b is small. This fact enables us to formulate all mechanical balance conditions 
in terms of a polygonal approximation of the cell array, thus allowing us to reduce the problem to a 
generalized "vertex model" [8j[9] . 

Remarkably, the forces given by ([T]) correspond to the mechanical energy in the form of the following 
simple Hamiltonian 

H({r l })=Y J H a [A ai {l ab }] (2) 

a 

where A a is the area of cell a, £ a b = \ r ij\ is the length of the interface between cells a and b and 
{(■ab} denotes the set of interfaces belonging to cell a. Both A a and £ a b's are defined in the polygonal 
approximation. This Hamiltonian is a generalization of the vertex models often used to describe epithelial 
sheet mechanics [8}[9j[26j. Pressure and tension are defined by considering the differential form of H: 

dH jn ^ dH , . 

(ab) a 

= ^T ah M ab + ^P a dA a 

(ab) a 

where we have define T ab = and P a = J^-. The (ab) sum runs over all edges, i.e. pairs of 

neighboring cells a, b. This tangent representation of mechanical energy expresses interfacial tension T 
and intracellular pressure P as conjugate variables to edge lengths and cell areas respectively. 
Mechanical equilibrium means that H is minimized with the respect to vertex positions 

*r-E*S~??ap-o <4) 

ah 

which defines the static force balance constraints. More generally, the dynamics of passive relaxation 
towards this equilibrium would be described by = — 8g gl ri ^ , where fj, is the "effective friction" 

constant. Our analysis will be based on the assumption that the cell layer is close to mechanical equi- 
librium in the sense of \Fi\ << meaning that the most of internal forces acting within the 
tissue are balanced. If there exists an unbalanced force that drives the physical motion, it is small in 
comparison to the forces that are balanced. 



The mechanical inverse problem. 

We can now inquire to what extent the knowledge that a given cell array geometry is in a mechanical 
equilibrium constrains the parameters P a , T ab describing the mechanical state of cells. We proceed by 
a simple count of mechanical constraints and of the free parameters for two cases i) a closed cell array, 
shown in Fig. [3]^ and ii) an open cell array, shown in Fig. [3j3. 
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Let us begin with the closed cell array and let v, e, c to be respectively the number of vertices, 
edges, and cells. For v vertices two dimensions, we have exactly n c = 2v — 3 mechanical constraints, 
where the extra three degrees of freedom are associated with global translation and rotation symmetries 
(alternatively, three constraints are redundant because the total force and total torque in the closed 
system are equal to zero). On the other hand, the number of unknown tension parameters is e, and the 
number of unknown pressures is c, so that the total number of parameters is n p = e + c. Our closed 
system, if we count exterior as an additional "cell", is topologically equivalent to a sphere so that Euler's 
theorem reads 

«-e + (c+l) = 2. (5) 

Combining this relation with the condition that vertices are points where three edges meet and each edge 
impinges on two vertices, that is 3i> = 2e, we obtain the result 

e + c = 2v + l. (6) 

This implies n p = n c + 4, which means that our unknown parameters can be determined up to four free 
constants. One of the latter is the arbitrary overall scale of T a f, and P a which cannot be constrained 
by the force balance conditions. (Note also since P a is only defined up to an additive constant, one can 
set the pressure in the exterior of the domain to zero.) Yet the good news is that the number of free 
constants is finite, while the number of nontrivial constraints scales with the number of cells! 

Repeating the counting procedure for the open system, one finds that e + c = 2v + b + 1, where b is 
the number of cells at the boundary of the domain. It follows that n p = n c + b + 1. Thus mechanical 
parameters are determined up to b + 1 free constants: we can still choose the overall scale while the 
additional b degrees of freedom may be regarded as the boundary conditions such as P a 's of the cells at 
the edge of the domain. Again, for a large array, because b ~ \fc while n p ~ c, the number of parameters 
and constraints is much larger than the number of free constants. 

To actually determine the T a b, P a parameters we use the fact that they appear only linearly in the 
force balance equations Q leading to a linear system for 

i> T = (T 1 ,...,T e ,P 1 ,...,P c ) (7) 

in the form 

MiP = C (8) 

with M being an n p x (n c + 1) matrix the 1st n c rows of which impose force balance conditions and the 
additional row imposing the scale, by constraining the average tension to be equal to one. Correspondingly 
the top n c entries of the column vector C are zero, while the bottom row C„ c +i = 1. 

The rectangular system ([8|) is solved by via pseudo-inverse 30 with the general solution of the form 
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^ = * + 1^ A ^ v ( 9 ) 



ith 



# = M _1 C, (10) 
M<f> v = (11) 

with M _1 being the pseudo-inverse of the rectangular matrix M and the amplitudes A v of the n z — 
n p — n c — 1 "zero modes" <f> v are the free parameters. 

Fixing the remaining n z degrees of freedom requires introducing additional constraints: e.g. one may 
have reasons to seek a solution which minimizes variation of P a 's or T a f,'s. In choosing such additional 
assumptions one may want to use all the information that one has for specific applications, as we shall 
do below. However, before proceeding to the applications we must consider the issue of error sensitivity. 
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Sensitivity of the inverse. 

Our approach to mechanical parameter inference is based on the observed geometry of the cell array. How 
sensitive are the results to the inaccuracy of vertex positions {r^}? Such an inaccuracy will inevitably 
arise in the process of imaging and image segmentation and even more importantly from the fact that cells 
fluctuate and our assumption that any particular configuration is in equilibrium, is at best approximate. 
To quantify the stability of the inverse we consider the effect of an arbitrary small perturbation in vertex 
positions, {Sfi}. Because the inhomogeneous term C in ^ is independent of cell geometry, the variation 
of parameters Sip in the response to positional error is given by 

M8il>+[^il>]6f=0 (12) 



Sip = LSf 



Or 



Sr. (13) 



Ideally the error response matrix L has small eigenvalues providing a relatively robust inverse. On 
the other hand, large eigenvalues of L would indicate high error sensitivity. These sensitive modes appear 
via the pseudoinverse matrix M^ 1 . A histogram of singular values of the matrix M~ l is shown in blue 
in Fig. [4] (for a closed system with « 800 cells) . One notes that a substantial fraction of modes have 
eigenvalues larger than one. As a result, small errors in positions can result in large error in inferred 
parameters. 

The simplest way to solve the sensitivity problem is to reduce the number of parameters. For example, 
as we shall argue below, in some contexts it may be reasonable to neglect variation in cell pressure and 
set P a = Pq which eliminates c parameters, reducing n p from 4c to 3c. In that case the mechanical 
constraint system given by ^ becomes overdetermined and can be solved only in the sense of least 
square minimization: i.e. minimization of 

Tr[{M'iP-C) T {M'4>-C)]. (14) 

The solution of the minimization problem is still given by the pseudo-inverse of the reduced rectangular 
matrix M' , the reduction being accomplished by eliminating constrained unknown parameters from in 
ip. Fig. [4] shows (in red) the distribution of singular values governing the sensitivity of the reduced or 
partial inverse problem. We note a substantial reduction in sensitivity. 

The partial inverse approach can tested in silico. To that end we consider a closed array of cells that 
appears in Fig. [4]^ and define cell geometry by minimizing elastic energy given by 

Hv({r}) = ^k ab (t ab -l) 2 . (15) 

(ab) 

with uniformly distributed k a b £ [0.7,1.3]. The absence of area terms imposes constant pressure. (The 
cell array is relaxed under toroidal boundary conditions to prevent a collapse into the zero tension ground 
state.) The vertex model parameters are computed via equation ([2]). These quantities are then compared 
to values obtained by applying the partial inverse algorithm to the vertex "data" {fi} corrupted by 
random noise {Sr^} with an r.m.s. variation of 5% of the average length of cell edge (see Fig. (|5|). The 
correlation coefficient between inferred and computed parameters is 0.852, which confirms the ability of 
our method to extract information from noisy data. 

We note that the "soft modes" which give rise to the sensitivity of the full inverse problem are 
quite interesting. The formulation of the minimally constrained problem is analogous to the isostatic 
systems studied in jamming transitions of amorphous solids [31]. These isostatic systems live on the 
boundary of Maxwell's criterion for rigidity, and much like amorphous solids, they must satisfy both 
the local and global rigidity conditions. In our mechanical inverse formulation, "rigidity" corresponds 
to a fully constrained set of mechanical (T a f, and P a ) parameters. Amusingly, local soft modes for the 
MI problem correspond to special local geometries: specifically, polygons that can be inscribed into 
circles (i.e. generalization of regular polygons) - a category which includes triangles of any shape. These 
interesting mathematical aspects of the problem will be discussed in a separate publication. 
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Results 

Mechanical differentiation of cells in the developing avian cochlea. 

During cochlear development, which takes place during the 1st two weeks of chick embryonic development, 
cells in the initially homogeneous two dimensional epithelial layer differentiate into pro-neural (hair-cell) 
and support cell fates [29]. The process is driven by the Delta/Notch-mediated cell-contact signaling [32] 
which causes lateral inhibition: cells which express Delta ligand on their surface prevent their immediate 
neighbors from doing the same. Expression of Delta is an early marker of the pro-neural fate of cells. Fig. 
1A presents an image of the cochlea epithelium, obtained by Goodyear and Richardson [29] at the stage 
of development shortly after the onset of differentiation. The two cell types already have a discernibly 
different morphology: pro-neural cells are somewhat smaller and have curved edges. This dimorphism is 
supported by direct labeling of specific pro-neural markers, shown in Fig. IB and demonstrated in [29] . 

Our goal is to infer, based on the analysis of the image in Fig. f A, the variation in the mechanical 
parameters between cells. The visible positive curvature associated with pro-neural cells suggests that 
they are under higher internal pressure. Can the Mechanical Inverse method determine pressure differ- 
entials between cells? Because our approach requires only positions of cellular vertices, it does not use 
the information provided by the interfacial curvatures which are readily measurable on the image. This 
additional information will be used as an a posteriori validation of the inferred results. 

To reduce the number of parameters we assume that interfacial tensions can be expressed as T a b = 
T a + Tb in terms of cortical tensions T a , Tb of adjacent cells. This reduces the number of parameters by 
e — c = 2c which is sufficient to render a robust partial inverse (in the sense of least squares), yielding 
T a and P a for every cell. Fig. [7] shows the distribution of inferred intracellular pressures and cortical 
tensions, for the two cell types. We see that pro-neural cells have on average higher tension and pressure. 
While pressure shows some correlation with cell area, there is no correlation between interfacial tension 
and its length. There, however is no reason to expect any specific correlation between these quantities. 
On the other hand, Laplace' Law predicts P a — Pb = n a bT a b which we are in a position to check directly, 
thanks to the fact that interfacial curvatures n a b are directly measurable on the images such Fig. [T]A Fig. 
[6]presents the "empirical" Laplace' Law obtained on the basis of the inferred P a — Pb and T a b. Because the 
Mechanical Inverse algorithm did not in any way use the interfacial curvature information, the fact that 
inferred parameters approximately obey the Laplace' Law provides a validation of the inverse method. 



Mechanical anisotropy at the onset of the ventral furrow formation in Drosophila. 

Ventral furrow formation in Drosophila is the first step of the gastrulation process and begins with the 
contraction of apical surfaces of cells along the ventral midline of the ellipsoidal monolayer of cells that 



comprise the embryo at that early stage of development 27 33 34 . Fig. 10A presents the ventral view 
of a Drosophila embryo at the beginning of this mechanical transformation. The high quality of these 
images makes it possible to attempt the Mechanical Inverse analysis. Since the process begins even before 
ccllularization is completed it is reasonable to assume that cells have the same internal pressure P a = P , 
allowing us to reduce the number of parameters enough to achieve a robust partial inverse and infer T a b 
for every cell boundary. We find a rather broad distribution of tensions (with the coefficient of variation 
« 0.3). 

Interestingly, comparing images separated by merely two minutes (Fig. [8| we found that the inferred 
T a b at the later time-slice exhibited statistically significant anisotropy with estimated tensions of cell 
interfaces along the AP axis being on average about 15% higher than those along the DV axis. The 
inferred increase in AP tension (relative to DV) is consistent with the laser ablation measurements made 
in the Wieschaus lab [7 27 . Yet, mechanical inverse inference gives information not only on the global, 
tissue- wide level, but also on the scale of a single cell and interface. The analysis also clearly demonstrates 
the ability to make specific predictions (for interfacial tensions) that can be directly tested by combining 
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high quality live imaging with UV pulsed laser ablation. 
Intercellular traction forces. 

The variation of tension from one interface to another implies the existence of traction forces acting 
between cells. This traction, or shear stress, must be entirely borne by the cadherins which bridge 
cellular membranes and connect actomyosin cortices of apposing cells [l] . In Fig. [9] we zoom in on an 
interface decomposing interfacial tension into the cortical tensions on the opposite sides of the interface 
T ao = T a +T = T a > + T D > , allowing for the possibility that the latter are not constant along the interface 
and vary as a function of position along the edge. (Here, T a and T a i refer to the two ends of the edge 
as shown in Fig. [9j) This transfer of tension from the cortical bundle in one cell to the other is possible 
because of cadherin mediated traction forces acting between cells. The total shear stress on the interface is 
Tab = (T a i — To) I L a b = {Tb — ) / L a b- In the Supplementary Information we show that because cortical 
tensions are constrained by the continuity conditions at cell "corners" they can be readily expressed 
in terms of interfacial tensions leading to the following simple expression for the traction force acting 
between cells a and b. 

1~ab=J (T ac - T ad + T bd ~ T bc ) (16) 

Fig. [l0) 3 shows inferred tractions calculated for the ventral furrow data. We observe a significant 
variability in tractions at different interfaces. Because traction forces stretch trans-cellular cadherin 
dimcrs, they may be physiologically important. Since at present there is no way of measuring them 
directly the possibility of indirect inference is particularly interesting. 



Discussion 

We have demonstrated that the readily visualized two dimensional network of cellular interfaces in an 
epithelial tissue holds, potentially, a wealth of information on the relative strength of mechanical stresses 
acting in the tissue. The main precondition is that the tissue is close to the mechanical equilibrium 
in which internal cytoskeletal forces are balanced by intercellular interactions. Any imbalance of forces 
corresponding to directed or fluctuating motion must be small by comparison with the balanced static 
component. Force balance is achieved by the suitable adjustment of cell geometries (parameterized by 
the positions of vertices) . Conversely we envision changes in tissue geometry to be driven adiabatically - 
i.e. without disruption of the mechanical equilibrium - by changes in cytoskeletal forces within cells. This 
picture is at once similar and dissimilar to the case of soap froths. The geometry of a soap froth 35 ■ 37 is 
also defined by the instantaneous force balance and changes adiabatically (when gas diffuses out of cells 
with higher internal pressure). Yet epithelial cells, in contrast to soap bubbles, can control interfacial 
tension by regulating myosin activity within actomyosin cortices and therefore can generate variation in 
tension on sub-cellular scale, even between different interfaces of the same cell. 



Our Mechanical Inverse method is fundamentally different from the Video Force Microscopy 28 . In 
contrast to our assumption that cytoskeletal forces are in an approximate instantaneous balance, VFM 
is based on the assumption that bulk forces acting within the tissue are balanced by viscosity: inverse 
is therefore based on the observed velocity of tissue motion. VFM employs finite element methods to 
define forces on a computational grid rather than underlying cells. The two methods are complementary 
in the sense that VFM provides information about the distribution of unbalanced bulk force which drives 
motion on the scale of the embryo, while our Mechanical Inverse is focused on the internal balance of 
forces in relation to cell geometry and its local changes. 

The proposed Mechanical Inverse method, converts clearly stated assumptions about the nature of 
cellular stresses into readily falsifiable predictions. Using the example of avian cochlea, we were able 
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to demonstrate that mechanical parameters inferred via the Mechanical Inverse satisfy non-trivial cross- 
checks provided by independent additional information (interfacial curvature measurements) read off the 
tissue images. Thus our approach is capable, in realistic applications, to infer mechanical parameters and 
to uncover interesting aspects of the internal state of the cell. By combining high quality live imaging 
with UV pulsed laser ablation, one will be able to put predictions for local interfacial tensions obtained 
via the Mechanical Inverse, to a rigorous experimental test. We note however, that the predictions do not 
have to be very accurate, to be useful. Even if inferred tensions each carry only one bit of information 
- i.e. identify interfaces with high or low tension - correlating tension with the observed level of myosin, 
cadherin and/or other proteins involved in regulation of cell mechanics could be extremely informative. 
(Since a large number of cells can be imaged and analyzed, the method is effectively "high throughput"!) 
In addition our approach allows to infer quantities such as inter-cellular traction forces (or shear stress), 
which may well be important for the stability of Adherens Junctions but cannot be directly measured 
by any presently available means. Hence we expect that further development, validation and application 
of the Mechanical Inverse method will leads to new insights into the molecular biology of epithelial cells 
and tissues. 
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Figure Legends 
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Figure 1. Avian cochlear epithelium (A) at the E9 stage of development just following differentiation 
of cells into hair-cell precursors and surrounding support cells (from Goodyear and Richardson [29]). 
Panel (B) shows the same tissue with stained (white) pro-neural cells. Panel (C) shows a computer 
generated segmentation of the raw image in (A) as a polygonal tiling which approximates cell geometry. 
The zoomed-in image (D) defines our parametrization of cell geometry in terms of vertex coordinates r, 
and of the mechanical state of the cell in terms of interfacial tensions T a b and hydrostatic pressures P a . 




Figure 2. Schematic representation of an edge between two cells in the tissue, comparing a curved 
interface (blue) with its approximation by a chord defining the edge in the polygonal representation of 
cells. Mechanical stress parameters are in red, and geometric quantities are labeled in black. Provided 
that the curvature of the interface K a b (and hence the angle 9) is small, the Young-Laplace equation 
P a — Pb — KabTab defines the force on the vertex i between cells a and b which obeys Eqn. (til). 
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Figure 3. Examples of computer generated closed (A) and open (B) cell arrays. Closed arrays provide 
an idealized context for denning the mechanical inverse problem, while the analysis of experimental 
data requires dealing with open arrays, corresponding to convex patches of cells defined by or within 
the field of view. 
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Figure 4. The distribution of singular values (which correspond to the square root of non-zero 
eigenvalues of [M T M] -1 ) for the M _1 matrix before (blue) and after (red) parameter reduction. Note 
that prior to parameter reduction there are a substantial fraction of eigenvalues > 10 which means that 
small errors in vertex positions are significantly amplified in solving the inverse problem. Large 
eigenvalues are effectively suppressed after parameter reduction. 




Figure 5. Scatter plot comparing actual T Q 6 describing the in silico cell array to the values inferred by 
the partial inverse algorithm applied to the vertex data corrupted by 5% random noise. The plot 
exhibits a correlation coefficient of 0.85 between the estimated and actual tensions. 
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Figure 6. Scatter plot comparing inferred pressure differential across an interface, P a — P&, with the 
product of inferred tension T a b and the measured curvature K a b of the same interface. Different colors 
distinguish results obtained from different images. The scatter plot exhibit a clear correlation between 
the two quantities, as expected from the Laplace' law P a — Pb = KabT a b- 
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Figure 7. Inferred tensions and pressures for the cochlear epithelium image shown in Fig. 1A. 
Hair-cell precursors and support cells correspond to red and blue dots respectively. Inferred pressure is 
plotted versus observed cell area and inferred tension versus edge length. Note systematically higher 
inferred pressure and tension in the hair-cells. 
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Figure 8. Histograms of inferred tension at the start of the ventral furrow formation. Red (blue) 
corresponds to cell edges at an angle above (below) 8 C — 7r/4 relative to the AP axis. Panels (A) and 
(B) correspond to respectively the 1st and the 3rd minutes of the furrow formation process. 
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Figure 9. Schematic decomposition of the effective interfacial tension into cortical tensions acting 
within cells: T a b = T a + Tb, Tb c = Tb + T c , etc. Because cytoskcletal cortexes of cells are crosslinked by 
cadherins via Adherence Junctions, indicated in blue, cortical stress can be transferred from one cell to 



another so that T a ^ T a > . The corresponding traction force (or shear stress) is given by Eqn. (16 1 
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Figure 10. Images of the ventral side of a Drosophila embryo 4min prior to ventral furrow 
invagination 27 . Panel (A) shows the polygonal tiling array defined by image segmentation. Panel 



(B) shows inferred tractions obtained from the partial inverse and Eq. ( |16| . Color indicates the 
magnitude of inferred traction with red (blue) being the relatively high (low) traction. The coefficient of 
variation of inferred traction is ~ 0.2. 



